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For the first time, we use the Gamma-ray Burst Monitor (GBM) on-board the Fermi satellite to 
search for sterile neutrino decay lines in the energy range 10-25 keV corresponding to sterile neutrino 
mass range 20-50 keV. This energy range has been out of reach of traditional X-ray satellites such 
as Chandra, Suzaku, XMM-Newton, and gamma-ray satellites such as INTEGRAL. Furthermore, 
the extremely wide field of view of the GBM opens a large fraction of the Milky Way dark matter 
halo to be probed. We start with 1601 days worth of GBM data, implement stringent data cuts, 
and perform two simple line search analyses on the reduced data; in the first, the line flux is limited 
without background modeling, and in the second, the background is modeled as a power-law. We 
find no significant excess lines in both our searches. We set new limits on sterile neutrino mixing 
angles, improving on previous limits by approximately an order of magnitude. Better understanding 
of detector and astrophysical backgrounds, as well as detector response, can further improve the 
limit. 


I. INTRODUCTION 

Right-handed neutral fermions (henceforth sterile 
neutrinos) arise in most implementations of the seesaw 
mechanism to generate neutrino masses, and yield ex¬ 
tremely rich phenomenology (for recent reviews, see, 
e-g-, Ul, Bl)- In particular, in the mass range of 1- 
100 keV, sterile neutrinos can be produced in suffi¬ 
cient quantities in the Early Universe to be a viable 
dark matter candidate. If the production proceeds 
via oscillations with active neutrinos, the momentum 
distribution of the sterile neutrino results in a warm 
dark matter candidate Q. However, other production 
mechanisms have been shown to result in sterile neu¬ 
trino dark matter that act very similarly to cold dark 
matter [3-Il|- Astrophysically, a dark matter sterile 
neutrino lies in the parameter space to be produced 
in core-collapse supernova cores providing a new 
mechanism for explosion 0 , as well as explaining the 
origin of strong neutron star kicks [HHIl- 

As a viable dark matter candidate, sterile neutrinos 
are stable on cosmological time scales, but neverthe¬ 
less they have decay channels that become interesting 
indirect detection targets for large concentrations of 
dark matter. The primary decay we target is the ra- 
diatively decay channel into an active neutrino and 
a photon. As the photon carries half of the sterile 
neutrino mass energy, the photon lies in the X-ray 
range. The signal is spectrally distinct from most ex¬ 
pected astrophysical backgrounds, since the signal line 
is broadened by the velocity dispersion of the dark 
matter particles. Coupled with the expected spatial 
morphology - spherical and centrally concentrated at 
the Galactic center - searches with X-rays can be a 
very powerful probe 0 - 


Many searches have been performed in the past, 
using X-ray satellites such as Chandra, Suzaku, and 
XMM-Newton\ observing a wide range of targets, from 
galaxy clusters, nearby galaxies, and dwarf satellites 
of the Milky Way, to the cosmic X-ray background. 
Most recently, an unexplained X-ray line was detected 
from a stack of ga laxy clusters as well as the An¬ 
dromeda galaxy (see also followup studies sup - 

porting and refuting these initial claims, e.g., [iTH^ l. 
which can be interpreted as the decay of 7 keV sterile 
neutrinos [0 . 

In this proceedings, we report initial results of us¬ 
ing the Gamma-ray Burst Monitor (GBM) onboard 
the Fermi satellite to search for X-ray lines arising 
from sterile neutrino decay. Among the advantages of 
using the GBM for this purpose include: (i) its all¬ 
sky coverage, which allows the entire Milky Way dark 
matter halo to be studied, and (ii) the energy range of 
the GBM, which fills a gap in energy coverage between 
X-ray satellites and gamma-ray satellites. We there¬ 
fore focus on the energy range = 10-25 keV cor¬ 
responding to sterile neutrino mass TOs = 20-50 keV, 
and explore the Milky Way because of its proximity 
and well-studied dark matter distribution. 


II. EXPECTED SIGNAL 

A. Intensity calcnlation 

The radiative decay of sterile neutrino has a decay 
rate of [3, [3| > 


/sin226l\ 


[ 10-10; 

V1 keV / 


eConf C141020.1 




2 


5th Fermi Symposium : Nagoya, Japan : 20-24 Oct, 2014 


where we have assumed a Majorana sterile neutrino. 

The photon intensity (number flux per solid angle) 
arising from sterile neutrino dark matter decay, from 
the direction of angle ip away from the Galactic Cen¬ 
ter, is 




dN 

dAdTdndE 


dTrmsTs 




( 2 ) 

dz dN \ 


where the first term in the brackets is the contribution 
from the Galactic halo and the second term is the con¬ 
tribution from extragalactic halos. Here, Tg = l/Tg is 
the sterile neutrino lifetime, pq = 0.4 GeV cm“^ is the 
local dark matter mass density, Rq = 8.5 kpc is the 
Sun’s distance to the galactic center, and dN/dE = 
S{E — msl2) is the photon spectrum. J{tp) is the so- 
called J-factor or boost factor, and is the integral of 
the dark matter mass density p in the Milky Way halo 
along the line-of-sight, 

JW = — dep{ilj,£), (3) 

PqRq Jo 

where imax is the outer limit of the dark matter halo. 
We assume the dark matter distribution is spherically 
symmetric about the Galactic Center, hence 

pi'ij,e) = p{rGc{i^,£)) (4) 

= pUR%-2tRQCOsil} + p\ . (5) 


We adopt £max = 250 kpc in this work. Although 
the value of £rnax differs depending on the adopted 
halo model, the contribution to from beyond 

^ 30 kpc are negligible. 

For the dark matter density profile p we adopt a 
NFW profile as our canonical profile. Although the 
Milky Way dark matter density profile at small radii 
remains uncertain, it is known well enough for ro¬ 
bust predictions of sterile neutrino decay signals on 
the scales of interest for the GBM. We adopt the fol¬ 
lowing generic dark matter halo profile, motivated by 
numerical simulations. 




/ r 

ri + (i?Q/i?,)“i 

\RqJ 

. l + {r/R,)°‘ _ 


(/3-7)/a 

1 

( 6 ) 


TABLE I: Dark matter profile parameters for widely 
adopted dark matter profiles in the literature. Our canon¬ 
ical profile is the NFW profile. 

Profile a P 'y Rs [pc] 

NFW 13 1 20 

cNFW 1 3 1.15 23.7 

Isothermal 2 2 0 3.5 



FIG. 1: The J-factor, Jijp), as a function of half open¬ 
ing angle ip relative to the Galactic Center, for four Milky 
Way dark matter halo profiles. Note the J-factor will be 
convolved by the GBM field of view, which is energy de¬ 
pendent but 40 degrees or more (see Figure[2]), which dra¬ 
matically reduces the difference between the dark matter 
profiles for the analysis. 


where parameters for commonly used profiles are sum¬ 
marized in Table m Another profile favored by recent 
simulations is the Einasto profile, 


p^”(r)=p0exp 


2 - Rq‘^ 


( 7 ) 


with as = 0.17 and scale radius Rs = 20 kpc. The 
difference in the J-factor between these profiles are 
shown in Figure [T] Here, the J-factors are shown as 
functions of the angle '0 away from the Galactic Cen¬ 
ter. When the J-factor are convolved by the GBM 
field of view, which is energy dependent but typically 
~ 40 degrees (see Figure [2]), the differences between 
dark matter profiles is dramatically reduced. 


B. Signal modeling for the GBM 

The Fermi-GBM consists of 12 Nal detectors and 2 
BGO detectors, the former covering energies 8 keV to 
1 MeV, and the latter covering 200 keV to 40 MeV. 
The Nal detectors are physically placed on the cor¬ 
ners and sides of the satellite. At any given time, 3 - 
4 Nal detectors view an Earth occultation, i.e., when 
the earth is within 60 degrees of the detector zenith. 
In the following, we perform a search using data from 
a single Nal detector, det-0. Det-0 is conveniently 
placed closest to the Fermi-LAT zenith (offset angle 
20.6°). Thus, as the LAT engages in survey mode, 
which are designed to maximize sky coverage, so does 
det-0. Analyses on additional detectors are forthcom¬ 
ing. 
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FIG. 2: The effective area for det-0 Nal detector versus 
photon arrival angle with respect to the detector normal. 
Three energy ranges utilized in the sterile neutrino search 
are shown. The occasional dips in the effective area are 
due to blockages from the satellite components. 


The response of det-0 are shown in Figures [5] and 
[21 where the effective area is shown as functions of 
the angle away from the detector normal, 6, and the 
photon energy. The angular dependence can be well- 
fitted with a Cosine function (shown by the solid line). 
The occasional dips are due to blockage from satellite 
components. 

One important consideration for the GBM is that 
it lacks photon-tracking capabilities, i.e., the photon 
count as a function of incident angle cannot be simply 
obtained. Earth occultation techniques can be effec¬ 
tively used to obtain directionality for point sources 
[^ . but remains difficult for diffuse sources. Fortu¬ 
nately, the lack of photon tracking ability is not very 
problematic for our sterile neutrino search, since the 
decay signal has a very large angular extent. However, 
it does mean it is difficult to accurately re-construct 
an intensity sky map (Eq. [5]) . We opt to simulate the 
instrument observable by properly modeling the sig¬ 
nal taking into account detector response. In this case, 
the instrumental observable is the photon count rate 
as a function of the Nal detector pointing direction. 

The photon counts is energy dependent and direc¬ 
tion dependent, i.e., Vi^j, where i labels the energy 
bin and j labels the detector pointing direction. The 
expected number of photons per observing time, Tj, 
from a particular detector pointing direction, is then 

^max 

^ = [ ' dE [ dm [ dE (8) 

^-^3 J J2tt J 

x\x{'4,,E)g(e,e) Heff(E,d)l, 



Energy [ keV ] 

FIG. 3: The same as Fig. [2]but plotted against energy, for 
three incident angles. The fits come from the angular fits 
for various energies. 

where Ef^^^ and E™™ are the boundaries of the i-th 
energy bin. We have integrated over the hemisphere 
the Nal detector points, i.e., over the detector zenith 
angle 9, and attribute all the photons to the pixel de¬ 
fined by pointing j. A position on the sky with an 
angle relative to the GC, ip, is thus related to the 
detector zenith angle and the pixel that the detector 
points at through ip —>• tp{9,j) . The pointing direc¬ 
tion of the detector is therefore defined by '0(0, j). The 
factor G{E, E) takes into account the energy resolu¬ 
tion of the Nal detector, which we model as a Gaus¬ 
sian with width given by the pre-launch calibrations 
[13, [il]- Lastly, Aea{E,9) is the Nal detector effec¬ 
tive area, which depends on energy and the detector 
zenith angle, as in Figures |2l|3] 

III. DATA REDUCTION 

We use a total of 1601 days worth of data, from 12- 
AUG-2008 to 31-DEC-2012. We use the CSPEC data 
with nominal 4.096 s time resolution and 128 channels 
in energy from 5 to 1402 keV. We then implement a 
wide range of data cuts to minimize a wide range of 
background contributions, which are summarized be¬ 
low. 

• LAT cuts: we use the Fermi-LAT cuts 
“LAT_CONFIG=l”, “LAT_MODE=5”, 

“DATA_QUAL=1”, “ROGK_ANGLE<50”, 

and “SAA=F”, where the first three conditions 
ensure the detector configurations and data 
qualities are suitable for scientific analysis, and 
the fourth condition ensures the Earth is not in 
the LAT’s field of view, which is approximately, 
although not exactly, the field of view of det-0 
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10-11 keV. Sin-2(ZTheta): 5.BK-11 



FIG. 4: All-sky counts maps in 10 - 11 keV energy range, showing the simulate dark matter map assuming the NFW 
profile (left) and the reduced data (right). The dark matter simulation assumes a sterile neutrino of mass 20 keV and 
mixing angle sin^ 26 = 5.8 x 10“^^. For det-0 only. 


(this is addressed later). The cuts also exclude 
epochs where the satellite is passing through 
the South Atlantic Anomaly (SAA), which 
has high cosmic-ray activity that significantly 
increases the radioactivity of the satellite. 

• Additional transient source cut: this removes 
epochs when the GBM detects transient sources, 
including gamma-ray bursts, direct cosmic-ray 
hits, solar flares, Galactic x-ray transients, and 
magnetospheric events. 

• Extended SAA cut: the LAT cuts do not com¬ 
pletely remove the effects of the SAA on the 
Nal detectors. The reason is that the satellite 
remains significantly radioactive even after lead¬ 
ing the SAA. We therefore remove the data col¬ 
lected in orbits that pass through the SAA suc¬ 
cessively. 

• Additional Earth cuts: we apply two addi¬ 
tional cuts to remove backgrounds related to the 
Earth. The first requires the angle between the 
Nal detector normal and the vector from the 
Earth center to be less than 50°, which reduces 
emissions from the Earth limb. The second re¬ 
quires the geomagnetic latitude to be less than 
|20|°, which reduces cosmic-ray induced vents. 

The reduced data product contains observed counts 
in 128 energy bins and 768 equal area sky pixels in 
healpix projection. The total live time after reduction 
equals ^ 4.6 x 10® seconds (^ 53 days). Despite the 
dramatic reduction in live time, the analysis is still 
systematic limited. In the energy range and region of 
interest, the total counts is more than ^ 10^ photons 
for each energy bin. 

Figure |4] shows the simulated dark matter photon 
count rate (left) and the reduced data count rate 
(right) for det-0. Both are shown with the same dy¬ 
namical range. Shown is the energy range 10-11 keV, 


which is our lower energy range. We observe a clear 
excess of photons towards the Galactic Genter direc¬ 
tion in the reduced data, which we interpret to be from 
astrophysical emissions from the Milky Way. Spec¬ 
tral analysis indicates the excess peaks at low ener¬ 
gies, and is dramatically reduced by 30-40 keV, when 
the data is dominated by cosmic-ray induced back¬ 
grounds. This is confirmed by the high-energy sky 
map showing small variations that trace the Earth 
magnetic field structure. 


IV. LINE SEARCH ANALYSIS 

Two line search strategies are implemented on the 
reduced data. The first is a conservative analysis 
based only on flux comparison, the second makes use 
of the spectral difference between signal (line) and 
background (dominated by a power-law within a small 
enough search energy window). 


A. Flux analysis 


The most conservative constraint on a sterile neu¬ 
trino decay amplitude is to require the decay sig¬ 
nal counts do not exceed the total measured photon 
counts. Thus, we do not make any assumptions of the 
detector background and astrophysical background. 

The comparison is made bin by bin, and consists of 
comparing the predicted signal 


Vi 


ROI 




dvij 

It;' 


to the observed photon counts, 


ROI 

di = ^ Nij , 
3 


(9) 


( 10 ) 
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FIG. 5: Constraints on the sterile neutrino mixing an¬ 
gle sin^(26) as a function of the sterile neutrino mass ms- 
Our new limits are shown by the dashed lines: the upper 
(black) line is from our conservative flux limit analysis, 
while the stronger lower (blue) line is from modeling the 
background as a power-law. The latter limit improves by 
preceding ones, obtained by observations of cosmic X-ray 
background (CXB) by HEAO-1 [s^, by about an order of 
magnitude. 


where Nij is the number of counts in energy bin i and 
pixel j measured by the det- 0 . 


B. Spectral analysis 


The flux analysis can be improved by modeling the 
background. We perform a spectral analysis which 
captures the different spectral shapes of the sterile 
neutrino decay signal and the backgrounds. 

The background is modeled as a power-law, 


dE 



M{E), 


( 11 ) 


7 o for all energy windows defined by Eq. (IT^ . We 
assume a Gaussian probability distribution function 
for each energy bin, giving the likelihood function as 


£(/3,7M) 


n 


'\/27rfJAeff 


(bj-did 


(13) 


where the product is taken over the energy bins, and 
CAeff is the systematic uncertainty of the effective 
area, 


CTAeff = 0.05(r'i -I- bi). (14) 


The present analysis is not limited by statistical un¬ 
certainty; rather, it is dominated by systematic un¬ 
certainties. We adopt a constant 5% uncertainty on 
the effective area as a conservative choice. The per 
degree of freedom for most energy windows is found to 
cluster between 1.1 and 1.4, although in some higher 
energy cases they can be as low as 0.5. We conclude 
that these findings justify the use of the local power- 
law assumption, as well as the 5% uncertainty in the 
effective area. 

Finally, the sterile neutrino decay line is added to 
the central energy bin for each energy window. The 
line signal has one free parameter, fs, the normal¬ 
ization, which scales linearly with the mixing angle 
8\t?29. Thus, the search has three free parameters, 
/3, 7 , and fs- The former two are treated as nuisance 
parameters, and fs is the parameter of interest. The 
likelihood is 


£{fs,K\d) 


n 


A/^TTfrAefF 


2“'Aeft^ 


(15) 


where k are the nuisance parameters. Finally, we use 
the profile likelihood analysis [ 2 ^ for treating nuisance 
parameters. In practice, this involves calculating the 
profile likelihood —ln£p(/s) for several fixed values of 
rus, where for each fs the —ln£ is minimized with 
respect to all other parameters k. 


V. STERILE NEUTRINO LIMITS 


where the normalization /3 and index 7 are left as free 
parameters. This assumption is made in a small en¬ 
ergy window defined by 

Max (iniin, i — Aui) < i < i -\- Aw, (12) 

using a fixed window of Aw = 5. This choice results 
in each energy width side being about 3 — 4 cr of the 
energy resolution. The energy window is truncated at 
*min = 6 , which corresponds to a central bin energy of 
9.4 keV. For line energies near the low energy cutoff, 
the energy window is thus asymmetric. 

We first perform a f^st to assess whether the 
assumption of a power-law is a good local background 
model. By way of minimizing the negative logarithm 
of the likelihood function, we find the best-fit do and 


We find no significant detections of lines in both 
our analyses. We determine the 95% C.L. one-sided 
upper limits on the amplitude, /g®, by requiring a 
2Aln£ = 2.71. The limits are shown in Figure [5] 
While the limits from the conservative flux analysis is 
weaker than previous limits, our spectral analysis limit 
results in an improvement of about a factor of ^ 10 
compared to those based on observations of cosmic X- 
ray background (CXB) by HEAO-1 [s^. The spectral 
analysis limit deteriorates at low energy, because the 
number of energy bins decreases, and also because the 
energy window is increasingly asymmetric. The latter 
in particular results in reduced ability to distinguish 
between the power-law background model and the line 
signal. 
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VI. CONCLUSIONS AND DISCUSSIONS 

We have used a GBM Nal detector (det-0) to set 
upper limits on sterile neutrino dark matter decays 
into mono-energetic photons. Two line analyses were 
implemented: the first is a conservative flux search, 
and the second is a spectral search assuming the back¬ 
ground can be modeled as a power-law over small en¬ 
ergy windows. In the energy range of 10-25 keV, cor¬ 
responding to sterile neutrino mass of 20-50 keV, our 
new upper limit is an improvement of about an order 
of magnitude compared to previous limits using the 
CXB with HEAO-1 data. 

The current analyses are dominated by systematic 
uncertainties primarily in the effective area. A better 
understanding of the GBM detector will therefore im¬ 


prove the limits presented in this proceeding. Work is 
currently underway to investigate other Nal detectors 
onboard the GBM. 
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